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Abstract 

Some dynamical properties for a dissipative time-dependent oval-shaped billiard are studied. The system is described in terms of 
a four-dimensional nonlinear mapping. Dissipation is introduced via inelastic collisions of the particle with the boundary, thus 
implying that the particle has a fractional loss of energy upon collision. The dissipation causes profound modifications in the 
dynamics of the particle as well as in the phase space of the non dissipative system. In particular, inelastic collisions can be 
assumed as an efficient mechanism to suppress Fermi acceleration of the particle. The dissipation also creates attractors in the 
system, including chaotic. We show that a slightly modification of the intensity of the damping coefficient yields a drastic and 
sudden destruction of the chaotic attractor, thus leading the system to experience a boundary crisis. We have characterized such 
a boundary crisis via a collision of the chaotic attractor with its own basin of attraction and confirmed that inelastic collisions do 
indeed suppress Fermi acceleration in two-dimensional time dependent billiards. 
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1. Introduction 

During the last decades many theoretical studies on dissipa- 
tive systems have been introduced in order to explain differ- 
ent physical phenomena in different fields of science including 
atomic and molecular physics^ 
namics 0,1111,01, optics 



X 



turbulent and fluid dy 
nanotechnology 111 ll Il2ll 



quantum and relativistic systems H 1 3L 114, 115H . Different pro- 
cedures have been used to describe such systems and two main 
different approaches are: (i) solving differential ordinary/partial 
equations or; (ii) using the so called billiard formalism. In prin- 
ciple, to chose procedure (i) or (ii) strongly depends on the type 
of system considered and possible existing symmetries. Case 

(i) are more likely devoted to systems where the external po- 
tential are smooth while case (ii) describes situations where the 
potential is null, say inside the boundary, and infinity outside 
of the boundary. The boundary identifies the position of this 
abrupt change. In this paper we shall concentrate to study case 

(ii) , i.e. a billiard system. 

A billiard consists of system in which one or many point-like 
particles move freely inside a closed region suffering specular 
reflections/collisions with the boundary. Billiards can be con- 
sidered one of the most attractive types of dynamical models in 
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the study of ergodic and mixed properties in Hamiltonian sys- 
tems 1 16]. From the mathematical point of view, a billiard is de- 
fined by a connected region Q c R D , with boundary dQ c R D ~ l 
which separates Q from its complement. If dQ = 8Q(f) the 
system has a time-dependent boundary and it can exchange en- 
ergy with the particle upon collision. Moreover, dissipation can 
be considered via different ways where the most common types 
used (i) drag force; (ii) damping coefficients. In the first case 
the particle loses energy/velocity as it were moving immersed 
in a fluid. For such a case, the dynamics is described by solv- 
ing differential equations lfl7ll . On the other hand, in the case 
(ii), the particle loses energy/velocity upon collision with the 
moving boundary. Thus the system is normally described using 
a billiard approach. It is know that depending on the combi- 
nation of initial conditions and control parameters, the phase 
space of such systems possess different structures. In the ab- 
sence of d issip ation, one kind of structure is the mixed type 
1 18 , Ijl 2(| 2ll 22, 23, 24, 25] where regular regions, such as in- 
variant tori and Kolmogorov-Arnold-Moser (KAM) islands are 
observed coexisting with chaotic seas. It is also well known in 
the literature that depending on the structure of the phase space 
the system can show or not a phenomenon called as Fermi ac- 
celeration, i.e., unlimited energy grow th j26Tl . Such a nomen- 
clature comes from Enrico Fermi [27] in 1949, as an attempt to 
explain the origin of cosmic ray acceleration. He proposed that 
such phenomenon was due to the interaction between charged 
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particles and time-dependent magnetic structures in the space. 
Since then the model has been modified and studied consider- 
ing different approaches. One of the most studied version of the 
problem is the Fermi-Ulam model (FUM) |28, 29]. Such model 
consists of a classical point-like particle moving between two 
rigid walls, one of them is assumed to be fixed and the other 
one moves according to a periodic function. In such system, 
Fermi acceleration is not observed since the phase space has a 
set of invariant tori limiting the size of the chaotic sea. How- 
ever, an alternative model was proposed by Pustylnikov l3(ill 
which is often called as bouncer model. Such system consists 
of a classical particle falling due to the action of a constant grav- 
itational field on a moving platform. One of the most important 
properties of this system is that depending on the combinations 
of both initial conditions and control parameters, the unlimited 
energy gain for a classical particle can be observed. It happens 
because there is no invariant tori limiting the size of the chaotic 
sea. 

When dissipation is taken into account, one can show that the 
mixed structure of the phase space present in the conservative 
case is destroyed. Then, an elliptic fixed point (generally sur- 
rounded by KAM islands) turns into a sink. Regions of chaotic 
seas might be replaced by chaotic attractors. Each one of these 
attractors has its own basin of attraction. Then, as a slight in- 
crease on the value of the damping coefficient, that is equivalent 
to reduce the power of dissipation, the chaotic attractor touch, 
even crosses, the line separating the basin of attraction of the 
chaotic attractor and the attracting fixed point (sink). Such be- 
haviour yields in a sudden destruction of the chaotic attractor. 
This destruction is called as a boundary crisis J31, 32, 33]. Af- 
ter the destruction, the chaotic attractor is replaced by a chaotic 
transient and its basin of attraction is destroyed, too. Addi- 
tionally, when dissipation is considered the behaviour of energy 
changes from unlimited to a constant plateau for long enough 
time. Thus, confirming that the phenomenon of Fermi acceler- 
ation is suppressed by dissipation 134 13511 . 

In the present paper we are interested in characterizing a 
boundary crisis in a time-dependent oval-shaped billiard. The 
paper is organized as follows. In Sec. |2]we describe all the need 
details to obtain the four-dimensional mapping that describe the 
dynamics of the system. Our numerical results are discussed in 
Sec. 12.11 Conclusion and acknowledgments are drawn in Sec. 

m 



2. The model and the mapping 

The two dimensional time-dependent oval billiard consists 
of a classical particle of mass m confined in and suffering col- 
lisions with a periodically moving boundary. The model is 2- 
dimensional (2-D) in the sense that it has two degrees of free- 
dom, however, the dimension of the phase space is defined 
as 2 X D. Based on this fact, we described the model us- 
ing a four dimensional and non linear map T(8„, a„,v„, t„) - 
(0„+i,a„ + i, V n +\,t„+{) where the dynamical variables are, re- 
spectively, the angular position of the particle; the angle that 
the trajectory of the particle does with the tangent line at the 
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Figure 1 : Plot of the boundary and a typical trajectory shown only for 6 snap- 
shots of the boundary. The control parameters used were e = 0.2, P = y = 1.0 
and r} = 0.05. 



position of the collision; the absolute velocity of the particle; 
and the instant of the hit with the boundary. The shape of 
the boundary is given in polar coordinates as Rb(0, e, rj, t) = 
[1 + ?7Cos(f)][l + ecos(2fJ)] where e is the circle's boundary de- 
formation and 77 is the amplitude of the time dependent pertur- 
bation and the sub-index b denotes boundary. Figure Q] shows a 
typical plot of the boundary and six collisions of a particle with 
the boundary. 

To construct the mapping, we start with an initial condition 
(9„,a„,V„,t n ). The Cartesian components of the boundary at 
the angular position (6„, t„) are 

X(0 n ,t n ) = [l+/7Cos(f„)][l + ecos(2fJ„)]cos(fJ„), (1) 
Y(6 n ,t n ) = [1 +77Cos(f„)][l +ecos(2(9„)]sin((9„). (2) 

The angle between the tangent of the boundary at the position 
(X(9 n ), Y(6 n )) measured with respect to the horizontal line is 



cf>„ = arctan 



Y'(0 n ,t n ) 



X'(9 n , ?„) 



(3) 



where the expressions for both X'(6„, t„) and Y'(8„, t„) are writ- 
ten as 



X'(0„,t„) = 
Y'(9„,t„) = 



dR(0 n ,t„) 

dR(0 n ,t„) 
d6 n 



cos(fJ„)-/?(fJ„,f„)sin(fJ„), (4) 
sm(0 n )+R(0„,t„)cos(0„), (5) 



with dR(9„, t„)/d8 n = -2e[l + ?7Cos(f„)] sin(26'„). Since the ex- 
pressions for <f> n and a„ are known, the angle of the trajectory 
of the particle measured with respect to the positive X-axis is 
(tp n + a„). Such information allows us to write the particle's 
velocity vector as 



~f„ = |V^|[cos(0„ + a„)i + sm((f>„ + a„)j] , 



(6) 



where i and j denote the unity vectors with respect to the X and 
Y axis, respectively. The position of the particle, as a function 



2 



of time, for t >t n , is given by 

X p (f) = X(8„,t„) + \P„\cos(<f>„ + a„)(t-t n ) , (7) 
Y p (t) = F(0„,f„) + |t„|sin(0„ + a„)(f-f„). (8) 

We stress the sub-index p denotes that such coordinates cor- 
respond to the particle. The distance of the particle mea- 
sured with respect to the origin of the coordinate system is 

given by R p (t) = ^Xj(t) + Yj(t) and 8 p at (X p (t), Y p (t)) is 
p = arctan[Yp(f)/^ r p(f)]- Therefore, the angular position at 
the next collision of the particle with the boundary, i.e. 9 n+ \, is 
numerically obtained by solving R p (8„+\,t n+l ) = Rb(8„ + ut„+i). 
It means that the position of the boundary is the same as the 
position of the particle at the instant of the collision. The time 
f„+i is obtained by evaluating the expression 



mvannat spanning curves 



tn+l — t n + 



VAX 2 + AY 2 



|t„| 



(9) 



where AX = X p (6 n+1 , t n+l ) - X{8„, t„) and AY = Y p {8 n+l , t n+l ) - 
Y(6 n ,t n ). To obtain the new velocity we should note that the 
referential frame of the boundary is moving. Since we are con- 
sidering inelastic collisions, the particle experiences a fractional 
loss of energy upon collision in both its normal and tangential 
components. Therefore, at the instant of collision, the following 
conditions must be matched 

V*„+i-7*„ +] = pV'n (10) 
V' n+x -t n+l = -yV' n -t n+l , (11) 

where the unitary tangent and normal vectors are 

7*„+i = cos(</>„ +I )/ + sin(0„ +I X/', (12) 
N n+ i = -sin(0„ + i)/ + cos(0„ + i)j. (13) 

/3 and y are damping coefficients, it means that the particle can 
loses velocity/energy upon collision in its normal component 
(y), tangential component (/3) or both. We consider both y e 
[0, 1] and /3 € [0, 1]. The completely inelastic collision happens 
when y = ft — and is not considered in this paper. On the other 
hand, when y — /3 = 1, corresponding to an elastic collision, all 
the results for the non-dissipative case are recovered. The upper 
prime indicates that the velocity of the particle is measured with 
respect to the moving boundary referential frame. At the new 
angular position 9 n+ \, we find that 



+ (i-#i& + i)^ w i 



(14) 



+ (l+y)f b (t n+ l)-% + l , (15) 

where ~^/,(f„+i) is the velocity of the boundary which is written 

as 

dRb(t n +i), 
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Figure 2: Phase space in the variables velocity and time for a time dependent 
oval billiard. The control parameters used were e = 0.2, yS = y = 1.0 and 
n = 0.05. We assume as fixed the initial conditions ao = tt/2 and 6q = 3n/2. 



with 



dRb{t n+ i) 



-7][l + ecos(2#„ +1 )] sin(f„ +1 ) 



dt n+ i 
Then we have 

|t„ +1 | = -y/o^+l • f n+I ) 2 + (f„ +1 -^„ +1 ) 2 
Finally, the angle a„+\ is written as 



a n+ \ = arctan 



(17) 



(18) 



(19) 



dt n+ i 



-[cos(0„ +1 )/ + sin(6>„ + iX/] 



(16) 



With this four dimensional mapping, we can explore now nu- 
merical results for the dynamics of the particle. 

2.1. Numerical Results 

In this section we discuss our numerical results. Just to re- 
mind, our main goal is to characterize a boundary crisis in a 
time-dependent oval-shaped billiard. To start, we show in Fig. 
[2] a typical phase space for a special set of initial conditions: 
ao = 7r/2 and 0q = 3^/2. For such combination of initial condi- 
tion and taken into account e = 0.2 and 77 = 0.05 the boundary 
has neutral curvature. With this particular choice of initial con- 
ditions, the phase space of the system is mixed. On the other 
hand, if we chose random ao and 80, the particle experiences 
the phenomenon of unlimited energy growth lf36Tl . 

We now consider the situation where both damping coeffi- 
cients p + 1 and y + \. We then keep fixed up to the end of 
the paper /3 = 0.25. It implies that there is a high dissipation 
along the tangent component of the particle's velocity. Results 
for different /? will be published elsewhere ||37m . The parameter 
y is considered from the order of y = 0.89. It is shown in Fig. 
Oa) the behavior of the attractors present in the system for the 
following combination of control parameters e - 0.2, (3 - 0.25, 
y - 0.8899 and 77 = 0.05. We can see a chaotic attractor and 
an attracting fixed point. Figure|3jb) shows their corresponding 
basin of attraction. The procedure used to construct the basin 




Figure 3: (Color online) (a) Attracting Fixed point indicated by a blue star 
(*) and chaotic attractor (set of black points); (b) Their corresponding basin of 
attraction. Black corresponds to the basin of attraction for the chaotic attractor 
and dark gray ( red) of the attracting fixed point. The control parameters used 
were e = 0.2, (i = 0.25, y = 0.8899 and n = 0.05. 

of attraction was divide both V and t into windows of 500 parts 
each, thus leading to a total of 2.5 x 10 s different initial con- 
ditions. Each initial condition was iterated up to n — 5 x 10 6 
collisions with the boundary. We see that only two attractors 
emerged for such combination of control parameters: sink and 
chaotic attractor. We stress that other attractors could in prin- 
ciple exist. If they exist however, their basin of attraction are 
too small to be obtained. It is clear that, after a very long num- 
ber of collisions of the particle with the boundary, the velocity 
of the particle does not grow unlimitedly. Consequently, no 
Fermi acceleration is observed and we conclude that introduc- 
tion of inelastic collisions worked out perfectly as a mechanism 
to suppress Fermi acceleration, as proposed in Ref. [38] for a 
stochastic 1-D system. 

Let us now go ahead with the characterization of the bound- 
ary crisis l3lLl32l 13311 . It is well known in the literature that a 
saddle fixed point, in the plane V x t has two kinds of manifolds: 
(a) stable and (b) unstable. The unstable manifolds are formed 
by a family of trajectories that turn away from the saddle fixed 
point. One of them can form the chaotic attractor (or visit the 
region of the chaotic attractor after the event of crisis), while 
the other one moves towards an attracting fixed point. These 
manifolds are obtained from the iteration of the map T with ap- 
propriate initial conditions. Similarly, the construction of stable 



Figure 4: (Color online) Basin of attraction for the chaotic attractor and at- 
tracting fixed point (sink). The region in black corresponds to the basin of 
attraction of the chaotic attractor; the region in dark gray ( red), denotes the 
basin of attraction of the attracting fixed point; light gray (green) in (a) iden- 
tifies the chaotic attractor. The control parameters used to construct the basin 
of attraction were e = 0.2, B = 0.25, r\ = 0.05 and y = 0.8899. The dissi- 
pation used in (a) the chaotic attractor were y = 0.8899 (before crisis, light 
gray (green)); (b) the chaotic transient were y = 0.8906 (after crisis, light gray 
(green)). 



manifolds are a little bit more complicated since the inverse of 
the mapping, say T l , must be obtained. The procedure for ob- 
taining the stable manifolds is the same as that one used for the 
unstable manifolds, however, instead of iterating the map T we 
must iterate its inverse T . Since the stable manifolds generate 
the border of the basin of attraction of the chaotic attractor and 
attracting sink, a boundary crisis happens when a chaotic attrac- 
tor touch the stable manifold due to a modification of the con- 
trol parameter. As a consequence, there is a sudden and drastic 
destruction of the chaotic attractor and its basin of attraction. 

It is shown in Fig. Ufa), two basins of attraction; one in black, 
corresponding to the basin of attraction of the chaotic attractor, 
and the other one in dark gray (red), denoting the basin of at- 
traction of the attracting fixed point, and the chaotic attractor 
marked by light gray (green). If we increase the value of the 
parameter y, which is equivalent to reduce the intensity of the 
dissipation, the two branches of the stable manifold touch, even 
crosses, the edges of the chaotic attractor, see Fig. H](b). Such 
behaviour is equivalent to a collision of the chaotic attractor 
with its own basin of attraction. Of course after the collision, 
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Figure 5: Behaviour of E X nfor different values ofy, as labeled in the figure. 
The control parameters used in the construction of the figure were p = 2, e = 
0.4, B=landrj = 0.001. 



there is a sudden destruction of chaotic attractor and its basin 
of attraction. We have observed such crisis for the control pa- 
rameters e = 0.2, p = 0.25, 77 = 0.05 and y = 0.8906. Other 
boundary crisis are observed for different combinations of con- 
trol parameters, too. After the boundary crisis, the entire plane 
Vxf(t mod (2n)) is the basin of attraction for a sink. There- 
fore, all initial conditions in such a region will converge to the 
sink. Additionally, for the regime of weak dissipation, the aver- 
age energy, £, = ^ 2" = o ^i,j> grows for short iteration number 
and suddenly it bends towards a regime of saturation for long 
enough values of n as can be seen in Fig. [5] Consequently, the 
mechanism of Fermi acceleration is suppressed in high as well 
as weak dissipation. 



3. Conclusion 

As a short remark, we have studied a classical version of a 
dissipative time-dependent oval-shaped billiard. The dissipa- 
tion was introduced via damping coefficients for both the nor- 
mal and tangential components of the particle's velocity. For 
the regime of high tangential dissipation, we characterized an 
event of boundary crisis. For the regime of weak dissipation, 
we have shown that the average energy remains constant for 
long enough time. Such result allows us to confirm that the in- 
troduction of inelastic collisions is sufficient to suppress Fermi 
acceleration since all the initial conditions will converge to at- 
tractors located at low velocity domain. 
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